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There are three basic mathematical tools we need, and then we can begin 
working on the physical implications of Schrodinger’s equation, which will 
take up the rest of the semester. So we start with a review of: Linear 
Algebra, Separation of Variables (SOV), and Probability. There will be 
no particular completeness to our discussions - at this stage, I want to 
emphasize those aspects of each of these that will prove useful. Pitfalls and 
caveats will be addressed as we encounter them in actual problems (where 
they won’t show up, so won’t be addressed). 

Try as one might, it is difficult to relate these three areas, so we will take 
each one in term, starting from the most familiar, and working to the least 
(that’s a matter of taste, so apologies in advance for the wrong ordering - in 
anticipation of our section on probability: How many different orderings of 
these three subjects are there? On average, then, how many people disagree 
with the current one?). 

We will begin with separation of variables, a simple technique that you have 
used recently - we will review the basic electrostatic interest in these solu¬ 
tions, and do a few examples of SOV applied to PDE’s other than Laplace’s 
equation (V 2 V = 0). 

1.1 Separation for the Laplace Problem 

Separation of variables refers to a class of techniques for probing solutions to 
partial differential equations (PDEs) by turning them into ordinary differen¬ 
tial equations (ODEs). We generally rely on some notion of uniqueness 1 for 

1 Uniqueness implying both a solution to the PDE and satisfaction of some boundary 
conditions. 
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our PDE - then the logic of SOV rests on the notion that finding a solution, 
even by this somewhat limited technique, suffices. 

There are two basic ffavors of SOV, additive and multiplicative - in either 
case, the style of argument is the same, so we’ll start with the algebraically 
simpler version. 


1.1.1 Additive Separation 


We take the Laplace problem: Find V ( x, y) satisfying V 2 V = 0 on some 
region, with some specified boundary conditions (for example, in E&M, we 
would have an implicit V —> 0 at spatial infinity condition, and maybe some 
conducting boundaries where we set V = 0). Working in two dimensions 
(for simplicity), we start with a separation ansatz- additive, in this case: 

V(x,y) = V x (x) + V y {y), (1.1) 


then running this through the Laplacian gives: 

d 2 14 d 2 V v 

-- H- y - = 0. 

dx 2 dy 2 


( 1 . 2 ) 


Now for the standard argument: V”(x) depends only on x just as Vy(y) 
depends only on y - they cannot be equal unless each is equal to a constant. 
Suppose we set V" = a, then we must have V” = —a. The solution to these 
ODEs is: 


V x (x) 


1 

2 


ax 2 + (3x + 5 


Vy(y) 


1 2 

~2 ay +72/ + P- 


(1.3) 


for arbitrary constants (/l, 5, 7 , p), to be used in imposing the boundary 
condition(s). The solution, for any a, can be written as: 


V(x,y ) 


- a (x 2 + y 2 ) + /3 x + 7 y + n 


(1.4) 


As far as electrostatics goes, this solution does not match our usual notion of 
potential - it doesn’t die at spatial infinity. That’s not always a deal-breaker, 
but it is true that additive separation is of limited utility in E&M (it plays a 
much larger role in, for example, Hamilton-Jacobi theory). Nevertheless, 
the pattern of: ansatz, argument, solution is the same as for the more 
immediately useful, and algebraically involved multiplicative separation. 
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1.1.2 Multiplicative Separation 


The difference here is in the starting ansatz - rather than additive single 
variable functions, we take the multiplicative combination: 


V(x,y) = V x (x ) V y (y), 


(1.5) 


then Laplace’s equation becomes 


d 2 V x 

dx 2 


Vy + V X 


d 2 Vy 

dy 2 


= 0 , 


( 1 . 6 ) 


and if we divide by V(x,y ) itself, we can write this in a form conducive to 
making the separation argument: 


V" V" 

^ + -dL = 0 . 
14 V,. 


(1.7) 


The first term depends only on x, the second only on y, so for a solution 
satisfying our functional assumptions, we must have both terms equal to a 
constant, call it a 2 (a convenient choice, as you will see in a moment) - for 
a complex. Then 


V" V" 

Y^L =a 2 = _ V _JL 

V x Vy 


( 1 . 8 ) 


will certainly have V 2 V = 0. The solutions here are familiar: 


V x = Ae ax + Be~ ax V y = Fe ioy + G e~ iay . (1.9) 


Now we see that for a generic complex number, both 14 and V y will be 
combinations of growing and decaying exponentials, as well as oscillatory 
sines and cosines. For example, if a is real, 14 is growing and decaying, 
while Vy is pure sine and cosine. For a = ia with a real, the roles are 
reversed. In general, we use the boundary conditions to choose a convenient 
expression for a. 


Example 

Consider a pair of infinite, grounded conducting sheets separated a distance 
d with a conductor connecting the two sheets held at 14 as shown in Fig¬ 
ure 1.1. What is the potential in between the plates? We know that the 
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potential satisfies Laplace’s equation in the region between the plates (no 
charge in there), and the boundary conditions are clear: 


V(x = 0, y) = V(x = d, y) = 0 V(x,y = 0) = Vq sin 


2 7 T X 


( 1 . 10 ) 


y 



Figure 1.1: Two parallel plates separated by a distance d and held at V = 0, 
the connecting plate is held at Vo sin(2 7 xx/d). 

We know the solution will be of the form (1.9), but how to choose a and the 
constants (A, B, F,G)? There is an additional, implicit boundary condition 
- we’d like the potential to go to zero in the “open” spatial direction, y —> 
oo - this tells us that we should set a = ia for a G 1R to get growing and 
decaying exponentials for the V y (y) solution: 

V x = Ae iax + Be~ iax V y = F e~ ay + Ge ay (1.11) 

and, furthermore, we should set G = 0 to eliminate the growing exponential. 
By introducing A = (A + B) and B = i (A — B), we can write V x {x ) in terms 
of sine and cosine. Our solution so far, then, is: 

V x = A cos(a x) + B sin(a x) V y = F e~ ay V = V x V y . (1.12) 

Now for the rest of the boundary conditions: take V x {x = 0) = 0 - that tells 
us that A = 0. Our solution is now much simpler, the full potential has the 
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form 

V = Q sin(ax) e~ ay (1.13) 

where we have introduced yet another rewrite of combinations of constants 
(Q = AF). The boundary condition V(x = d,y) = 0 is interesting, this 
tells us that sin(ad) = 0, not a statement about the obvious constant Q, 
but rather, a constraint on the allowed wavelengths: 

sin(acf) = 0 —> ad = hit (1.14) 

for integer n. Then the only a’s that satisfy the boundary condition are 
a = still an infinite family, but labelled by the integer index n. So, we 
have 

V n (x,y) = Q n sm(mrx/d)e~ n7ry/d , (1.15) 

and we put the subscript n here to remind us that any integer n provides 
a solution. This type of boundary condition, one that limits the allowed 
“wavelength”, is what leads to the mathematical representation of quanti¬ 
zation when applied to the separable solutions of Schrodinger’s equation. 

We have one boundary left - because Laplace’s equation is a linear PDE, 
sums of solutions are still solutions, and we can make a general solution 
out of our n-indexed ones. Each term satisfies Laplace’s equation and the 
boundary conditions at x = 0, d and y —> oo, so the sum does as well: 

QO _ 

V(x,y) = ^Q n sin(^-^) e ^. (1.16) 

n =1 

As for the final boundary condition: V(x,y = 0) = Vq sin(2 n x/d), we can 
take Q n = 0 for all n except for n = 2. We know that the solution is unique, 
so this is certainly valid. But it begs the question, what about a less well- 
adapted driving potential? That is a question for next time, as it is tied to 
the decomposition of functions in sin and cosine “bases”, a linear algebra 
issue. 

The final answer: 

V(x,y) = Vo s™(1.17) 

This solution is just the potential at the y = 0 plate, decaying into the 
region between the plates - equal contours are shown in Figure 1.2. 
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Figure 1.2: The potential between the plates for the setup shown in Fig¬ 
ure 1.1. 


There is not much more to say about separation of variables - the best 
approach is to do as many examples as you can get your hands on - that 
makes some of the choices that appear unmotivated at first (why should we 
make the integration constant a 2 ? Why do we pick the y solution to be 
growing and decaying?) more reasonable. 


1.2 Separation for Other PDEs 

We will be applying the SOV technique to Schrodinger’s equation, and that 
is not going to be much different than the Laplacian example from E&M. 
The approach is independent of PDE, although the PDE must be linear (not 
strictly speaking true, but the superposition principle we invoked in writ¬ 
ing (1.16) relied on linearity) - just to show some of the places where SOV 
comes up, we can look at the Helmholtz equation and the wave equation. 
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1.2.1 Helmholtz Equation 

Suppose we consider the same problem as above, but for the Helmholtz 
PDE: 

(V 2 + n 2 ) V = 0 

V(x,y ->• oo) —> 0, 
(L18) 

we are imagining a slight modification to electrostatics (actually amounts to 
a theory of electrostatics with a massive photon, but forget about that for 
now). 

First we’ll find the ODE solutions to the Helmholtz PDE via V(x,y) = 

V x {x)V y {y ): 

V" V" 

t + i + » =0 - (L19) 

We see immediately that we will again get oscillatory and exponential solu¬ 
tions, since we want the solution in the increasing y direction to decay, we’ll 

y/f 

call y = a 2 , with a real. With the constant /i 2 in place, we don’t need to 
set the x term equal to —a 2 , we have additional freedom. We know that the 

yff 9 

boundary conditions for x are easily imposed for sines, so set -y- = — /3 for 
/? real. Then the constraint we must impose to solve the Helmholtz PDE 
is: — /3 2 + a 2 + /r 2 = 0. Imposing the boundary conditions at x = 0, d gives 
back (3 = ^ as before, so that 

V x = A sin(^p) . (1.20) 

Now, for the V y equation, we are given /j, and we have (3 2 = 
the growing and decaying exponentials have the form: 

y y — p e V n2 TT 2 /d 2 -fj, 2 y _|_ Q e -\J n2 ^/(P-iFy' 

There are interesting elements to this solution, most notably, 
for which V y is actually oscillatory (n < yr ) - ignoring that for 
if we assume the above terms for V y are indeed growing and 
n = 2, we have the solution: 

V(x, y) = Rosin e ~yf^W^y^ (1.22) 

which just has a different fundamental decay length when compared to the 
Laplace form for electrostatics. 


n 2 7r 2 /d 2 , so 

( 1 . 21 ) 

a cutoff in n 
the moment, 
decaying for 


V(x, y = 0, d) = 0 V(x,y = 0) = Vq sin 


2 7 T X 
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1.2.2 The Wave Equation 


Separation of variables is not limited to purely spatial problems - we can 
mix in time. For the wave equation in one spatial dimension (call it x), 
defined in a medium with fundamental speed v, we have 

1 d 2 u(x,t) d 2 u(x,t ) 

v 2 dt 2 dx 2 


Take a multiplicative ansatz: u(x, t ) = u x (x) ut(t), then by the now familiar 
procedure of inputting and dividing by u: 



// 

t. 


u t 



V*X 


(1.24) 


and again, setting both terms equal to a constant: 



// 

t 


Ut 


= a 


< 

Ux 


= —a 


(1.25) 


What is interesting here is that both the spatial solution u x and the temporal 
ut can be oscillatory - that is not true for the Laplace equation separable 
solutions we saw above. Without worrying about boundary conditions, we 
have the easy solution: 


ut = A cos (avt) u x = B cos(ax) 


(1.26) 


so we can write u(x, t) = uq cos(a x) cos(a vt) - of course, just as obvious is 
the solution u(x,t) = uq sin(ax) sin (avt). Since the wave equation is lin¬ 
ear, the sum of these two solutions is also a solution, and it is an instructive 
one: 

u(x, t ) + u(x, t ) = uq (cos(a x) cos(a vt) + uq sin(a x) sin(a v t)) 

= uq cos(cc (x — vt)), 

precisely a plane wave. Of course, in general, we would be given spatial 
boundary conditions and/or an initial waveform at t = 0. 


(1.27) 


1.2.3 Poisson’s Equation 

The utility of separation is not limited to source-free equations - consider 
Poisson’s equation for electrostatic potential in the presence of source charge 
density p: 


(1.28) 
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We must be given a p that is itself appropriately separable. The simplest 
possible such distribution would be a constant, but we can imagine more 
interesting charge distributions. As an example, suppose we are in spherical 
coordinates, and we have a spherically symmetric charge density: p(r, 6, q i) = 
p(r), depending only on our distance from the origin. Poisson’s equation 
becomes, in these variables 

/ 1 8 ( 2 dV\ 1 d ( . <9R\ 1 <9 2 R\ p(r) 

\r 2 dr \ dr ) r 2 sin# dd \ Sm d9 ) r 2 sin 2 6 dcj) 2 ) eo 

(1.29) 

The separation ansatz now becomes V(r,6,(j)) = V r (r) Vq(9) V^(^>), and we 
can see what will happen - to satisfy the functional dependence assumption, 
we will have constant Vg and V^, leaving us with just the radial ODE: 


1 <1 
r 2 dr 



P(r) 

eo 


(1.30) 


Suppose we take p{r) = po a constant, then the solution to the above is: 


t / / \ . P Por 2 . . 

W(r) = aH--—. (1.31) 

r 6 e 0 

Notice that the first two terms are just solutions to V 2 V = 0 (an overall 
constant and the potential outside a spherically symmetric distribution of 
charge) - if we are inside a uniformly charged sphere, with r = 0 included 
in the domain, then we must set (3 = 0. The constant cc, of course, can be 
set to zero, and we are left with the usual potential for a uniformly charged 
sphere. 
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Homework 

Reading: Griffiths “Introduction to Electrodynamics”, pp. 127-136. 
Problem 1.1 

Solve Laplace's equation V 2 V = 0 using separation: V(x, y) = V x (x) V y {y) 
for the following two-dimensional boundary conditions: V(0,y) = 

V(d, y) = 0, V(x,d) = 0 and V(x,0) = Vq sin(^p). 



Physically, we are describing the electrostatic potential inside a square that 
has grounded lines on three sides, and a spatially-varying potential on the 
fourth side. 


Problem 1.2 

The heat equation in one (spatial) dimension reads 

1 du(x,t) _ d 2 u{x,t.) n QrA 

k dt 8x 2 ’ 1 j 

where u(x,t ) is the temperature of an object with “thermal conductivity” k 
(a material property). For a complete solution, we must specify two bound- 
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ary conditions (since this equation is second order in x ) and an “initial” 
condition (the value of u(x,t = 0), for example). 

a. Using separation of variables: u(x,t ) = u x (x)ut(t), and separation 
constant —a 2 , write down the separated heat equation and solve for u x and 

u t . 

b. Suppose we have a material that goes from x = 0tox = d-we 
set the temperature on both ends of the material to zero - i.e. u(0,t) = 
u(d,t) = 0. Initially, at time t = 0, we have temperature: u(x,0) = 
«o **(”)■ Using your separation solution, find the particular solution for 
this set of boundary and initial conditions. 
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